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ABSTRACT 

Quantum technology relies on the utilization of resources, like quantum coherence and entanglement, which allow quantum 
information and computation processing. This achievement is however jeopardized by the detrimental effects of the environment 
surrounding any quantum system, so that finding strategies to protect quantum resources is essential. Non-Markovian and 
structured environments are useful tools to this aim. Here we show how a simple environmental architecture made of two 
coupled lossy cavities enables a switch between Markovian and non-Markovian regimes for the dynamics of a qubit embedded 
in one of the cavity. Furthermore, qubit coherence can be indefinitely preserved if the cavity without qubit is perfect. We 
then focus on entanglement control of two independent qubits locally subject to such an engineered environment and discuss 
its feasibility in the framework of circuit quantum electrodynamics. With up-to-date experimental parameters, we show that 
our architecture allows entanglement lifetimes orders of magnitude longer than the spentaneous lifetime without local cavity 
couplings. This cavity-based architecture is straightforwardly extendable to many qubits for scalability. 


Introduction 

Entangled states are not only an existing natural form of com¬ 
pound systems in the quantum world, but also a basic resource 
for quantum information technology.*"^ Due to the unavoid¬ 
able coupling of a quantum system to the surrounding envi¬ 
ronment, quantum entanglement is subject to decay and can 
even vanish abruptly, a phenomenon known as early-stage dis¬ 
entanglement or entanglement sudden death.Harnessing 
entanglement dynamics and preventing entanglement from 
disappearing until the time a quantum task can be completed 
is thus a key challenge towards the feasibility of reliable quan¬ 
tum processing. '"*’ 

So far, a lot of researches have been devoted to entangle¬ 
ment manipulation and protection. A pure maximally en¬ 
tangled state can be obtained from decohered (partially en¬ 
tangled mixed) states*®"^" provided that there exist a large 
number of identically decohered states, which however will 
not work if the entanglement amount in these states is small. 
In situations where several particles are coupled to a common 
environment and the governing Hamiltonian is highly sym¬ 
metric, there may appear a decoherence-free subspace that 
does not evolve in time:^'"^^ however, in this decoherence- 
free subspace only a certain kind of entangled state can be 
decoupled from the influence of the environment.^"*’^^ The 
quantum Zeno effect^® can also be employed to manipulate 
decoherence process but, to prevent considerable degradation 
of entanglement, special measurements should be performed 
very frequently at equal time intervals.^"*’By encoding each 
physical qubit of a manyqubit system onto a logical one com¬ 


prising several physical qubits,an appropriate reversal 
procedure can be applied to correct the error induced by de¬ 
coherence after a multiqubit measurement that learns what 
error possibly occurred. Yet, as has been shown,^* in some 
cases this method can indeed delay entanglement degrada¬ 
tion but in other cases it leads to sudden disentanglement for 
states that otherwise disentangle only asymptotically. The 
possibility to preserve entanglement via dynamical decou¬ 
pling pulse sequences has been also theoretically investigated 
recently for finite-dimensional or harmonic quantum envi- 
ronments^^"^^ and for solid state quantum systems suffering 
random telegraph or 1 // noise,but these procedures can 
be demanding from a practical point of view. 

In general, environments with memory (so-called non- 
Markovian) suitably strucmred constitute a useful tool for 
protecting quantum superpositions and therefore the entan¬ 
glement of composite systems.^’It is nowadays well- 
known that independent qubits locally interacting with their 
non-Markovian environments can exhibit revivals of entan¬ 
glement, both spontaneously during the dynamics^^’and 
on-demand by local operations.These revivals, albeit pro¬ 
longing the utilization time of entanglement, however even¬ 
tually decay. In several situations, the energy dissipations of 
individual subsystems of a composite system are responsible 
for disentanglement. Therefore, methods that can trap system 
excited-state population would be effective for entanglement 
preservation. A stationary entanglement of two independent 
atoms can be in principle achieved in photonic crystals or 
photonic-band-gap materials"*^’if they are structured so as 


to inhibit spontaneous emission of individual atoms. This 
spontaneous emission suppression induced by a photonic crys¬ 
tal has been so far verihed experimentally for a single quantum 
dot"*^ and its practical utilization for a multi-qubit assembly 
appears far from being reached. Quantum interference can 
also be exploited to quench spontaneous emission in atomic 
systems^**’and hence used to protect two-atom entangle¬ 
ment provided that three levels of the atoms can be used.^^ 
Since the energy dissipations originate from excited state 
component of an entangled state, a reduction of the weight 
of excited-state by prior weak measurement on the system 
before interacting with the environment followed by a reversal 
measurement after the time-evolution proves to be an efficient 
strategy to enhance the entanglement.However, the suc¬ 
cess of this measurement-based strategy is always conditional 
(probability less than one).^^“^^ It was shown that steady-state 
entanglement can be generated if two qubits share a com¬ 
mon environment,^"*’^® interact each other®^ and are far from 
thermal equilibrium.®^"®^ It has been also demonstrated that 
non-Markovianity may support the formation of stationary en¬ 
tanglement in a non-dissipative pure dephasing environment 
provided that the subsystems are mutually coupled.®® 

Separated, independent two-level quantum systems at ther¬ 
mal equilibrium, locally interacting with their own environ¬ 
ments, are however the preferable elements of a quantum 
hardware in order to accomplish the individual control re¬ 
quired for quantum information processing.**’ *® Therefore, 
proposals of strategies to strongly shield quantum resources 
from decay are essential within such a conhguration. Here we 
address this issue by looking for an environmental architecture 
as simple as possible which is able to achieve this aim and 
at the same time realizable by current experimental technolo¬ 
gies. In particular, we consider a qubit embedded in a cavity 
which is in turn coupled to a second cavity and show that this 
basic structure is able to enable transitions from Markovian to 
non-Markovian regimes for the dynamics of the qubit just by 
adjusting the coupling between the two cavities. Remarkably, 
under suitable initial conditions, this engineered environment 
is able to efficiently preserve qubit coherence and, when ex¬ 
tended to the case of two noninteracting separated qubits, 
quantum entanglement. We finally discuss the effectiveness 
of our cavity-based architecture by considering experimental 
parameters typical of circuit quantum electrodynamics,*®’®* 
where this scheme can hnd its natural implementation. 

Results 

Our analysis is divided into two parts. The first one is ded¬ 
icated to the single-qubit architecture which shall permit us 
to investigate the dynamics of quantum coherence and its 
sensitivity to decay. The second part treats the two-qubit ar¬ 
chitecture for exploring to which extent the time of existence 
of quantum entanglement can be prolonged with respect to its 
natural disappearance time without the proposed engineered 
environment. 



Figure 1. Scheme of the single-qubit architecture. A 

two-level atom (qubit) is embedded in a cavity Ci which is in 
turn coupled to a second cavity C 2 by a coupling strength J. 
Both cavities are taken at zero temperature and can lose 
photons. 

Single-qubit coherence preservation 

The global system is made of a two-level atom (qubit) inside 
a lossy cavity Ci which in turn interacts with another cavity 
C 2 , as depicted in Fig. 1. The Hamiltonian of the qubit and 
two cavities is given hy {h = 1) 

M = (Ct)o/2) C7^-f Oi -f Cp2^2^2 “b O’—-f C7+) 

-|-7(d|d2+ 0102 )’ 

where = |1) (1| — |0) (0| is a Pauli operator for the qubit 
with transition frequency cuo, d± are the raising and lowering 
operators of the qubit, di (d|) and d 2 (dj) *he annihilation 
(creation) operators of cavities Ci and C 2 which sustain modes 
with frequency coi and 0 ) 2 , respectively. The parameter K 
denotes the coupling of the qubit with cavity Ci and J the 
coupling between the two cavities. We take coi = (02 = CO and, 
in order to consider both resonant and non-resonant qubit- 
Ci interactions, Oo = ft) -b 5 with 5 being the qubit-cavity 
detuning. Taking the dissipations of the two cavities into 
account, the density operator p (f) of the atom plus the cavities 
obeys the following master equation®® 

Pit) = -i[H,p(f)] 

2 

- E y [ai^a«P(f)-2a„p(f)a*-bp(f)aJa„](2) 

n=l 

where p(f) = dp{t)/dt and F] (r 2 ) denotes the photon decay 
rate of cavity Ci (C 2 ). The rate r„/2 physically represents 
the bandwidth of the Lorentzian frequency spectral density of 
the cavity C„, which is not a perfect single-mode cavity.®® A 
cavity with a high quality factor will have a narrow bandwidth 
and therefore a small photon decay rate. Weak and strong 
coupling regimes for the qubit-Ci interaction can be then 
individuated by the conditions K<T\/A and tf> F i/4.**’®® 
Let us suppose the qubit is initially in the excited state 11) 
and both cavities in the vacuum states |00), so that the overall 
initial state is p(0) = |100) (100|, where the hrst, second and 
third element correspond to the qubit, cavity Ci and cavity C 2 , 
respectively. Since there exist at most one excitation in the 
total system at any time of evolution, we can make the ansatz 
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for p(f) in the form 

Pit) = il-Xit))\wit)) (Wf)| +A(f) 1000) (0001, (3) 

where 0 < A(f) < 1 with A( 0 ) =0 and |v^(f)) = hit) | 100 ) + 
Cl (f) | 010 ) +C 2 (f) | 001 ) with hiO) = 1 and ci ( 0 ) = C 2 ( 0 ) = 0 . 
It is convenient to introduce the unnormalized state vector^®’ 

\Yit)) ^ yr^xwiwo) 

= h(f)|100)+Q(f)|010)+?2(f)|001), (4) 

where hit) = •^/r^^X(t)/r(f) represents the probability ampli¬ 
tude of the qubit and c„(f) = -^/l — A(f)c„(f) (n = 1 , 2 ) that 
of the cavities being in their excited states. In terms of this 
unnormalized state vector we then get 


p (f) = I v/(f)) (wit) I + ^ (f) 1000) (0001. ( 5 ) 

The time-dependent amplitudes hit), ci(f), C 2 (f) of Eq. (4) 
are determined by a set of differential equations as 

(a)-l-5)/!(f)-|-K'ci(f), 

^cu- Ci(f)-|-K-/!(f)-|-7c2(f), 

^CU- ^r2^ C2(f)+-/Ci(f). (6) 

The above differential equations can be solved by means of 
standard Laplace transformations combined with numerical 
simulations to obtain the reduced density operators of the 
atom as well as of each of the cavities. In particular, in the 
basis { 11 ), | 0 )} the density matrix evolution of the qubit can 
be cast as 


,dhit) 

dt 

.dci it) 
dt 

.dciit) 

dt 



MrPll(O) 

2;poi(o) 


ztpioio) \ 

1 -M,pii( 0 ) J’ 


(7) 


where Ut and Zt are functions of the time t (see Methods). 

An intuitive quantification of quantum coherence is based to 
the off-diagonal elements of the desired quantum state, being 
these related to the basic property of quantum interference. 
Indeed, it has been recently shown®^ that the functional 

nt)= L \PiM, ( 8 ) 


where p,; (f) (i 7 ^ j) are the off-diagonal elements of the sys¬ 
tem density matrix, satisfies the physical requirements which 
make it a proper coherence measure.®^ In the following, we 
adopt ^(f) as quantifier of the qubit coherence and explore 
how to preserve and even trap it under various conditions. To 
this aim, we first consider the resonant atom-cavity interaction 
and then discuss the effects of detuning on the dynamics of 
coherence. 

Suppose the qubit is initially prepared in the state |<)>(0)) = 
a|0)-fj3|l) (with |a|2-f-|j3|2=i), namely, -^( 0 ) = 2|aj3|, 
then at time f > 0 the coherence becomes “^(f) = 2|aj3/!(f)|. 


Focusing on the weak coupling between the qubit and the 
cavity Ci with K = 0.24ri, we plot the dynamics of coher¬ 
ence in Fig. 2(a). In this case, the qubit exhibits a Markovian 
dynamics with an asymptotical decay of the coherence in the 
absence of the cavity C 2 (with J = 0). However, by introduc¬ 
ing the cavity C 2 with a sufficiently large coupling strength, 
quantum coherence undergoes non-Markovian dynamics with 
oscillations. Moreover, it is readily observed that the decay 
of coherence can be greatly inhibited by increasing the Ci- 
C 2 coupling strength J. On the other hand, if the coupling 
between the atom and the cavity Ci is initially in the strong 
regime with the occurrence of coherence collapses and re¬ 
vivals, the increasing of the C 1 -C 2 coupling strength J can 
drive the non-Markovian dynamics of the qubit to the Marko¬ 
vian one and then back to the non-Markovian one, as shown 
in Fig. 2(b). This behavior is individuated by the suppres¬ 
sion and the successive reactivation of oscillations during the 
dynamics. It is worth noting that, although the qubit can expe¬ 
rience non-Markovian dynamics again for large enough J, the 
non-Markovian dynamics curve is different from the original 
one for 7 = 0 in the sense that the oscillations arise before 
the coherence decays to zero. In general, the coupling of 
C 1 -C 2 can enhance the quantum coherence also in the strong 
coupling regime between the qubit and the cavity C\. 

The oscillations of coherence, in clear contrast to the mono¬ 
tonic smooth decay in the Markovian regime, constitute a 
sufficient condition to signify the presence of memory effects 
in the system dynamics, being due to information backflow 
from the environment to the quantum system.®^ The degree 
of a non-Markovian process, the so-called non-Markovianity, 
can be quantified by different suitable measures. ®®“22 \Yg 
adopt here the non-Markovianity measure which exploits the 
dynamics of the trace distance between two initially different 
states Pi ( 0 ) and P 2 ( 0 ) of an open system to assess their distin- 
guishability.®® A Markovian evolution can never increase the 
trace distance, hence nonmonotonicity of the latter would im¬ 
ply a non-Markovian character of the system dynamics. Based 
on this concept, the non-Markovianity can be quantified by a 
measure jV defined as®® 

max f a[t,pii0),p2i0)]dt, (9) 

Pl(0),p2{0)7a>0 

where ( 7 [t,pi( 0 ),p 2 ( 0 )] = £/Z)[pi(f),p 2 (f)]/£/f is the rate 
of change of the trace distance, which is defin ed as 
D[pi(f),p 2 (f)] = (l/2)Tr|pi(t)-p 2 (f)|, with |A| = y/X^'X. 
By virtue of jY, we plot in Fig. 3 the non-Markovianity of 
the qubit dynamics for the conditions considered in Fig. 2(a)- 
(b). We see that if the qubit is initially weakly coupled to 
the cavity Ci (k = 0 . 24 r 1) its dynamics can undergo a tran¬ 
sition from Markovian (M^ = 0) to non-Markovian (jY > 0) 
regimes by increasing the coupling strengths J between the 
two cavities. On the other hand, for strong qubit-cavity cou¬ 
pling (k = O.TFi), the non-Markovian dynamics occurring for 
7 = 0 turns into Markovian and then back to non-Markovian 
by increasing 7. We mention that such a behavior has been 
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Figure 2. Coherence '^(f) of the qubit as a function of the scaled time Fif for different coupling strengths J between the two 
cavities for (a) K = 0.24ri, r 2 = O.SFi and (b) K = 0.4Fi, F 2 = 0.5Fi. The qubit is initially prepared in the state |0(O)) with 
a = p = 1 /v/2 and resonant with the cavity (detuning 5 = 0). The plots in panels (c) and (d) display the coherence trapping 
for a perfect cavity (F 2 = 0) with K = 0.24F 1 and K — 0.4Fi, respectively. 



Figure 3. Non-Markovianity quantifier ^ of equation (9) 
of the qubit dynamics as a function of y/Fi for weak 
{k — 0.24Fi) and strong {k — 0.4Fi) coupling regimes to 
cavity Ci and a fixed decay rate F 2 = 0.5F 1 of the cavity C 2 . 


also observed in a different structured system where a qubit 
simultaneously interacts with two coupled lossy cavities. 

Trapping qubit coherence in the long-time limit is a use¬ 
ful dynamical feature for itself that shall play a role for the 
preservation of quantum entanglement to be treated in the 
next section. We indeed find that the use of coupled cavities 
can achieve this result if the cavity C 2 is perfect, that is F 2 = 0 
(no photon leakage). The plots in Figure 2(c)-(d) demonstrate 
the coherence trapping in the long-time limit for both weak 
and strong coupling regimes between the qubit and the cavity 



Figure 4. Density plots of coherence “^(f) of the qubit as 
functions of detuning 5 and the scaled time Fif for (a) 

K = 0.24Fi, F 2 = 0.2Fi, y = 0.5Fi; (b) K = 0.24Fi, 

F 2 = 0.2Fi, y = Fi; (c) ff = 0.4Fi, F 2 = 0.5Fi, J = 0.5Fi; 
(d) K = 0.4Fi, F 2 = 0.5Fi, y = Fi. The initial state of the 
qubit is maximally entangled (a ~ /3 = 1 /v/2)- The values 
of the coherence are within the range: 0 1. 
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Figure 5. Coherence ‘^(f) of the qubit as a function of the scaled time Fif for different values of the detuning 5 in the case 
when the cavity C 2 is perfect, that is r 2 = 0. The qubit-Ci and the C 1 -C 2 coupling strengths are, respectively, (a) K = 0.24ri, 
J = O.SFi; (b) K = 0.4ri, J = O.SFi. Out of resonance (5 > 0) no coherence trapping is achievable. 


Cl for different coupling strengths J between the two cavities. 
This behavior can be explained by noticing that there exists 
a bound (decoherence-free) state of the qubit and the cavity 
C 2 of the form |!/_) = J \ 10) — K |01), with J and K being the 
C 1 -C 2 and qubit-Ci coupling strengths. Being this state free 
from decay, once the reduced initial state of the qubit and 
the cavity C 2 contains a nonzero component of this bound 
state I v/_), a long-living quantum coherence for the qubit can 
be obtained. For the initial state |O(0)) = a |000) -f p \ 100) 
of the qubit and two cavities here considered and r 2 = 0, 
the coherence defined in Eq. (8) gets the asymptotic value 
CXD^ = 2\apJ^/{J^ + ff^)|, which increases with J for 
a given K. We further point out that the cavity Ci acts as a 
catalyst of the entanglement for the hybrid qubit-C 2 system, 
in perfect analogy to the stationary entanglement exhibited by 
two qubits embedded in a common cavity.In the latter case, 
in fact, the cavity mediates the interaction between the two 
qubits and performs as an entanglement catalyst for them. 

We now discuss the effect of non-resonant qubit-Ci interac¬ 
tion (5 ^ 0) on the dynamics of coherence. In Figure 4(a)-(d), 
we display the density plots of the coherence as functions of 
detuning 5 = (Oq — CO and rescaled time Ff for both weak and 
strong couplings. One observes that when 5 departures from 
zero, the decay of coherence speeds up achieving the fastest 
decay around 5 =J. It is interesting to highlight the role of the 
cavity-cavity coupling parameter 7 as a benchmark for having 
the fastest decay during the dynamics under the non-resonant 
condition. For larger detuning tending to the dispersive regime 
(5 ^ k), the decay of coherence is instead strongly slowed 
down.^^ However, as shown in Fig. 5, stationary coherence 
is forbidden out of resonance when the cavity Ci is perfect. 
Since our main aim is the long-time preservation of quantum 
coherence and thus of entanglement, in the following we only 
focus on the condition of resonance between qubit and cavity 
frequencies. 


Two-qubit entanglement preservation 

So far, we have studied the manipulation of coherence dynam¬ 
ics of a qubit via an adjustment of coupling strength between 
two cavities. We now extend this architecture to explore the 
possibility to harness and preserve the entanglement of two 
independent qubits, labeled as A and B. We thus consider A 
(B) interacts locally with cavity Cia (Cib) which is in turn 
coupled to cavity C 24 (€ 23 ) with coupling strength (Jb), as 
illustrated in Fig. 6 . That is, we have two independent dynam¬ 
ics with each one consisting of a qubit j {j = A,B) and two 
coupled cavities C\j-C 2 j- The total Hamiltonian is then given 
by the sum of the two independent Hamiltonians, namely, 
H = HjHj, where each Hj is the single-qubit Hamiltonian of 
Eq. (1). Denoting with F ij (F 2 ;) the decay rate of cavity Cij 
iC 2 j), we shall assume F^ = Fib = F as the unit of the other 
parameters. 

As known for the case of independent subsystems, the 
complete dynamics of the two-qubit system can be obtained 
by knowing that of each qubit interacting with its own en- 
vironment."*^*’^^ By means of the single-qubit evolution, 
we can construct the evolved density matrix of the two 
atoms, whose elements in the standard computational basis 
{|l) = |ll),|2)^|10),|3)^|01),|4) = |00)}are 

Pii(f) = ufufpn{0) 

P22(f) = uf{l-uf)pu{0) + i4P22{0) 

P33{t) = {^-uf)ufpil{0) + ufp22,{Q) 

P44(f) = (l-M)')(l-Mf)pil(0)-f (1 -m)^)P22(0) 

+ (1 ~ “f )P33(0) -|-P44(0) 

Pl4(f) = P4lit)=4^fPl4iO) 

P23(0 = Phit)=4zf*p23i0), (10) 

where pim{0) are the density matrix elements of the two-qubit 
initial state and m/,z/ are the time-dependent functions of Eq. 
(7). 
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Figure 6. Scheme of the two-quhit architecture. Two 

independent qubits A and B, initially entangled, are locally 
embedded in a cavity Cij which is in turn coupled to a second 
cavity C 2 j by a coupling strength Jj (j = A,B). 


We consider the qubits initially in an entangled state of the 
form Iv^(O)) = a |00) +j3 |11) (|ap + |j3p = 1). As is known, 
this type of entangled states with |)3| > |a| suffers from en¬ 
tanglement sudden death when each atom locally interacts 
with a dissipative environment.^"® As far as non-Markovian 
environments are concerned, partial revivals of entanglement 
can occur^^’^^^^’^^^^ typically after asymptotically decaying 
to zero or after a finite dark period of complete disappear¬ 
ance. It would be useful in practical applications that the 
non-Markovian oscillations can occur when the entanglement 
still retain a relatively large value. With our cavity-based 
architecture, on the one hand we show that the Markovian 
dynamics of entanglement in the weak coupling regime be¬ 
tween the atoms and the corresponding cavities (i.e., Cia and 
Cib) can be turned into non-Markovian one by increasing 
the coupling strengths between the cavities C 1 A-C 24 and (or) 
C\B-C 2 B', on the other hand, we find that the appearance of 
entanglement revivals can be shifted to earlier times. We 
employ the concurrence^^ to quantify the entanglement (see 
Methods), which for the two-qubit evolved state of Eq. (10) 
reads ‘^ab( 0 = 2max{0, |pi 4 (f)l “ A/P 22 (f)P 33 (f)}- Notice 
that the concurrence of the Bell-like initial state |i/(0)) is 
%b( 0) = 2|a)3|. In Fig. 7(a) we plot the dynamics of concur¬ 
rence %B(f) in the weak coupling regime between the two 
qubits with their corresponding cavities with Ka = Kb — 0.2r 
(Eia = Ei^ = r has been assumed). For two-qubit initial states 
with a — Y^l/ 10 , j5 — ^^9/10, the entanglement experiences 
sudden death without coupled cavities (Ja =Jb = 0). By incor¬ 
porating the additional cavities with relatively small coupling 
strength, e.g., Ja = O.SF and Jb — F, the concurrence still 
undergoes a Markovian decay but the time of entanglement 
disappearance is prolonged. Increasing the coupling strengths 
Ja, Jb of the relevant cavities drives the entanglement dynam¬ 
ics from Markovian regime to non-Markovian one. Moreover, 



Figure 7. The dynamics of concurrence for different 
coupling strengths Ja and Jb in (a) weak qubit-cavity 
coupling regimes with Ka = Kb = 0.2r and (b) strong 
qubit-cavity coupling regimes with Ka — Kb — 2r. The initial 
state weights are chosen as (a) a = -^ 1 / 10 , j3 = -^9/10 and 
(b) a = -^1/3, p = yJ'JjT), while in both cases 
r 2 A = ^ 2 B = 0.2r. The inset in (b) shows the long-time 
dynamics of concurrence for Ja = TT and Jb = ST. Panels (c) 
and (d) show the density plots of the two-qubit concurrence 
as a function of J (Ja = Jb = J is here assumed) and scaled 
time Ft, the others parameters being as in panels (a) and (b), 
respectively. The values of the concurrence in the density 
plots range within the interval; (c) 0 0 . 6 ; (d) 


the entanglement revivals after decay happen shortly after the 
evolution when the entanglement still has a large value. In 
general, the concurrences are enhanced pronouncedly with 
Ja and Jb- A comprehensive picture of the dynamics of con¬ 
currence as a function of coupling strength J is shown in Fig. 
7(c) where we have assumed Ja = Jb = J- In Fig. 7(b) we 
plot the dynamics of ’^Asit) in the strong coupling regime 
between qubit j and its cavity Cij with Ka = Kb = 2 r for 
which the two-qubit dynamics is already non-Markovian in 
absence of cavity coupling, namely the entanglement can re¬ 
vive after dark periods. Remarkably, the figure shows that 
when the coupling Jj between Cij and C 2 ; is activated and 
gradually increased in each location, multiple transitions from 
non-Markovian to Markovian dynamics surface. We point out 
that the entanglement dynamics within the non-Markovian 
regime exhibit different qualitative behaviors with respect to 
the first time when entanglement oscillates. For instance, for 
Ja=Jb = 3r, the non-Markovian entanglement oscillations 
(revivals) happen after its disappearance, while when Ja = 4r 
and Jb — SF the entanglement oscillates before its sudden 
death. These dynamical features are clearly displayed in Fig. 
7(d). 

As expected according to the results obtained before on 
the single-qubit coherence, a steady concurrence arises in the 
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Figure 8. The dynamics of concurrence for different coupling strengths and Tg in the presence of ideal coupled cavities 
C 2 A and C 2 B with r 2 A = r 2 g = 0 for (a) Ka = 0.2r, Kb = 0.3r and (b) Ka = Kb= 2r. The other parameters are chosen as 
a = 1 / 1 / 3 , j3 = •/2/3. The inset in (b) shows the short time dynamics of concurrence. 


long-time limit if the secondary cavities Cja, Qb do not lose 
photons, i.e., r2A = r2B = 0 . Fig. 8(a) shows the dynam¬ 
ics of concurrence for qubits coupled to their cavities with 
strengths Ka = 0 . 2 r, Kb — 0 . 3 r. We can readily see that, in 
absence of coupling with the secondary cavities (Ja =Jb — 0 ), 
the entanglement disappear at a finite time without any re¬ 
vival. Contrarily, if the local couplings C\j-C 2 j are switched 
on and increased, the entanglement does not vanish at a finite 
time any more and reaches a steady value after undergoing 
non-Markovian oscillations. Furthermore, the steady value 
of concurrence is proportional to the local cavity coupling 
strengths Ja, Jb- In Fig- 8(b), the concurrence dynamics for 
Ka = Kb = 2r is plotted under which the two-qubit entan¬ 
glement experiences non-Markovian features, that is revivals 
after dark periods, already in absence of coupled cavities, as 
shown by the black solid curve for Ja = Jb — 0 . Of course, 
in this case the entanglement eventually decays to zero. On 
the contrary, by adjusting suitable nonzero values of the local 
cavity couplings a considerable amount of entanglement can 
be trapped. As a peculiar qualitative dynamical feature, we 
highlight that the entanglement can revive and then be frozen 
after a finite dark period time of complete disappearance (e.g., 
see the inset of Fig. 8(b), for the short-time dynamics with 
Ja = 2 r, Jb — 3 r and also Ja — Jb — 3 r). We finally point 
out that the the amount of preserved entanglement depends 
on the choice of the initial state (i.e., on the initial amount 
of entanglement) of the two qubits. As displayed in Fig. 9, 
the less initial entanglement, the less entanglement is in gen¬ 
eral maintained in the ideal case of r2A = r2B = 0 . However, 
since there is not a direct proportionality between the evolved 
concurrence %iB{t) and its initial value %b( 0), the maximal 
values of concurrence do not exactly appear at a = 1 / 
(corresponding to maximal initial entanglement) at any time 
in the evolution, as instead one could expect. It can be then 
observed that nonzero entanglement trapping is achieved for 


a > 0 . 2 . 

Experimental paramaters 

We conclude our study by discussing the experimental feasibil¬ 
ity of the cavity-based architecture here proposed for the two- 
qubit assembly. Due to its cavity quantum electrodynamics 
characteristics, our engineered environment finds its natural re¬ 
alization in the well-established framework of circuit quantum 
electrodynamics (cQED) with transmon qubits and coplanar 
waveguide cavities.®^’The entangled qubits can be ini¬ 
tialized by using the standard technique of a transmission-line 
resonator as a quantum bus.®'*’®*' Initial Bell-like states as 
the one we have considered here can be currently prepared 
with very high fidelity.®" Considering up-to-date experimental 
parameters'®"^® applied to our global system of Fig. 6, the 
average photon decay rate for the cavity C\j {j = A,B) con¬ 
taining the qubit is Fiy S [1 MHz, 10 MHz], while the average 
photon lifetime for the high quality factor cavity € 2 ] is T 2 « 55 
which implies r 2 ; « lO^^MHz S [lO^^ri^, 10"*'ri;]. 
The qubit-cavity interaction intensity Kj and the cavity-cavity 
coupling strength Jj are usually of the same order of mag¬ 
nitude, with typical values Kj ^ Jj S [1 MHz, 100 MHz] = 
[0. IF ij, lOFij]. The typical cavity frequency is O) ~ 2;r x 10 
GHz®'* while the qubit transition frequency can be arbitrarily 
adjusted in order to be resonant with the cavity frequency. 
The above experimental parameters put our system under the 
condition Kj <C CO which guarantees the validity of the rotating 
wave approximation (RWA) for the qubit-cavity interaction 
here considered in the Hamiltonian of equation (1). 

In order to assess the extent of entanglement preservation 
expected under these experimental conditions, we can ana¬ 
lyze the concurrence evolution under the same parameters 
of Fig. 8(a) for Kj, Jj, which are already within the exper¬ 
imental values, but with r 2 A = r 2 B = r 2 = lO^^F, lO^^F 
instead of being zero (ideal case), where F = F^ = Fib € 
[1 MHz, 10 MHz]. The natural estimated disappearance time 
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Discussion 
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Figure 9. The concurrence as a function of the two-qubit 
initial state parameter a and the scaled time Tf for Ka = 0.2r, 
Kb = 0.3r, Ja = O.Sr, Jb = T and T 2 A = ^ 2 b = 0. The 
parameter a quantifies the initial entanglement according to 
the concurrence ^ab (0) = 2|aj3| = 2|a|-^/l — |ap. 


of entanglement in absence of coupling between the cavities 
(Ja =JB = 0)ist — 6.69/r G [669 ns, 6.69 /is], as seen from 
Fig. 8(a). When considering the experimental achievable de¬ 
cay rates for the cavities Cij, we find that the entanglement is 
expected to be preserved until times t* orders of magnitude 
longer than t, as shown in Table 1 . In the case of higher quality 
factors for the cavities C 2 ;, such that the photon decay rate is 
of the order of r 2 = lO^^T, the entanglement can last even 
until the order of the seconds. These results provide a clear 
evidence of the practical powerful of our simple two-qubit 
architecture in significantly extending quantum entanglement 
lifetime for the implementation of given entanglement-based 
quantum tasks and algorithms. 

It is worth to mention that nowadays cQED technologies 
are also able to create a qubit-cavity coupling strength com¬ 
parable to the cavity frequency, thus entering the so-called 
ultra-strong coupling regime.®^ In that case the RWA is to 
be relaxed and the counter-rotating terms in the qubit-cavity 
interaction have to be taken into account. According to known 
results for the single qubit evolution beyond the RWA,^^ it 
appears that the main effect of the counter-rotating terms in 
the Rabi Hamiltonian is the photon creation from vacuum 
under dephasing noise, which in turns induces a bit-flip error 
in the qubit evolution. This photon creation would be instead 
suppressed in the presence of dissipative (damping) mech¬ 
anisms.®'^ Since our cavity-based architecture is subject to 
amplitude damping noise, the qualitative long-time dynamics 
of quantum coherence and thus of entanglement are expected 
not to be significantly modified with respect to the case when 
RWA is retained. These argumentations stimulate a detailed 
study of the performance of our proposed architecture under 
the ultra-strong coupling regime out of RWA, to be addressed 
elsewhere. 


In this work, we have analyzed the possibility to manipulate 
and maintain quantum coherence and entanglement of quan¬ 
tum systems by means of a simple yet effective cavity-based 
engineered environment. In particular, we have seen how an 
environmental architecture made of two coupled lossy cavi¬ 
ties enables a switch between Markovian and non-Markovian 
regimes for the dynamics of a qubit (artificial atom) embedded 
in one of the cavity. This feature possesses an intrinsic interest 
in the context of controlling memory effects of open quantum 
systems. Moreover, if the cavity without qubit has a small 
photon leakage with respect to the other one, qubit coherence 
can be efficiently maintained. 

We mention that our cavity-based architecture for the single 
qubit can be viewed as the physical realization of a photonic 
band gap for the qubit,®^ inhibiting its spontaneous emission. 
This property, then extended to the case of two independent 
qubits locally subject to such an engineered environment, 
has allowed us to show that quantum entanglement can be 
robustly shielded from decay, reaching a steady-state entan¬ 
glement in the limit of perfect cavities. The emergence of this 
steady-state entanglement within our proposed architecture 
confirms the mechanism of entanglement preservation when 
the qubit-environment interaction is dissipative: namely, the 
simultaneous existence of a bound state between the qubit 
and its local environment and of a non-Markovian dynamics 
for the qubit.^® We remark that this condition is here shown 
to be efficiently approximated within current experimental 
parameters such as to maintain a substantial fraction of the 
entanglement initially shared between the qubits during the 
evolution. Moreover, we highlight that this goal is achieved 
even if the local reservoir (cavity) embedding the qubit is 
memoryless, thanks to the exploitation of an additional good- 
quality cavity suitably coupled to the first one. Specifically, 
we have found that, by suitably adjusting the control parame¬ 
ter constituted by this local cavity coupling, the entanglement 
between the separated qubits can be exploited for times orders 
of magnitude longer than the natural time of its disappearance 
in absence of the cavity coupling. These times are expected 
to be long enough to perform various quantum tasks. 

Our long-living quantum entanglement scheme, besides 
its simplicity, is straightforwardly extendable to many qubits, 
thus fulfilling the scalability requirement for complex quan¬ 
tum information and computation protocols. The fact that the 
qubits are independent and noninteracting also allows for the 
desirable individual operations on each constituent of a quan¬ 
tum hardware. The results of this work provide new insights 
regarding the control of the fundamental non-Markovian char¬ 
acter of open quantum system dynamics and pave the way to 
further experimental developments towards the realization of 
devices able to preserve quantum resources. 
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Tz/F 

Ta/F = Jb/Y = 0.5 

Ja/Y = 0.5, Jb/Y = 1 

10-^ 

f* =454/rG [45.4/ts,454/ts] 

f* =974/rG [97.4/IS,974/ts] 

O 

1 

u; 

t* = 4481/F e [448 /ts,4.48 ms] 

t* = 9686/F e [0.967 ms,9.67 ms] 


Table 1 . Estimates of the experimental entanglement lifetimes t* for different values of the second cavities decay rates r 2 and 
the local cavity couplings Ja, Jb- These values are to be compared with the natural entanglement lifetime without cavity 
coupling, t S [669 ns,6.69 /is]. The reference unit T G [1 MHz, 10 MHz]. 


Methods 

Functions of the single qubit density matrix 

Let us denote with the inverse Laplace trans¬ 

form of L{s). Then, the functions Ut and zt appearing in Eq. 
(7) are expressed as 

Mr = krP, Zt=^^'^{F{s)/G{s)}{t), 

where 

F{s) = -AJ^-{2s + 2i(0 + Yx)(2s + 2m + T 2 ){ll) 
G{s) = 2K'^(2s-|-2/a)-l-r2)-f [s-l-i(5-l-a))] 

X A~ A~ -1-2(5' -1- /cu)r2 
-fT I (2s -f 2i(0 -|- r 2 )}. 

Entanglement quantification by concurrence 

Entanglement for an arbitrary state Pab of two qubits is quan¬ 
tified by concurrence^’ 

%B = ‘i^(PAB)=max{0,Vzr-V^-Vl3-v^}> (12) 

where (/ = 1,..., 4) are the eigenvalues in decreasing order 
of the matrix pAB(o’y 0 Gy)p^g(^ay 0 a,,), with a,, denoting the 
second Pauli matrix and p^g corresponding to the complex 
conjugate of the two-qubit density matrix Pab in the canonical 
computational basis {[11), [10), jOl) , [00)}. 
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